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' Summary. Theory of expectation values is presented as an alternative to S-matrix 

, theory for quantum fields. This change of emphasis is conditioned by a transition 

^SJ ' from the accelerator physics to astrophysics and cosmology. The issues discussed are 

the time-loop formalism, the Schwinger-Keldysh diagrams, the effective action, the 
vacuum currents, and the effect of particle creation. 

Oh' 

D ' Introduction 

High-energy physics will probably have to undergo major changes. The accel- 
erators will cease being its experimental base, and it will become a part of 
^ I astrophysics. Simultaneously, the S-matrix will cease being the central object 

of high-energy theory because the emphasis on this object is entirely owing 
to the accelerator setting of the problem. If there is a background radiation 

■ that originates from some initial state in the past, then where is the S-matrix 
here? Astrophysics and cosmology offer the evolution problems rather than 
the scattering problems. The gravitational collapse is a typical initial-value 

I problem. It is such by its physical setting irrespective of whether the state of 

■ the system is classical or quantum. The nature of measurement also changes. 
No final state is prepared. One measures observables like temperatures or 
mechanical deflections and subjects these measurements to a statistical treat- 
ment to obtain the value of the observable. This means that one measures 

■ expectation values in the given initial state. S-matrix theory should give way 
to expectation-value theory. 

There is a proof that accelerator physics is dead: Gabriele Veneziano is 
leaving CERN for College de France. At this historic moment, my mission is 
to convert him into a new faith. The present preaching consists of 4 lectures: 

1. Formal aspects of expectation- value theory. 

2. The in-vacuum state and Schwinger-Keldysh diagrams. 



* The course of 4 lectures given at College de France in May 2006. 
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3. The effective action. 

4. Vacuum currents and the effect of particle creation. 

Literature to Lectures 1 and 2 is in [T]-[Tn]. Additional literature to Lecture 3 
is in [17]-[1I] and to Lecture 4 in [42]-[56l. 

1 Formal Aspects of Expectation- Value Theory 
Vocabulary 

In these lectures, 

0' (1.1) 

denotes the quantum field. It is an operator function on a given differentiable 
manifold (referred to below as the base manifold), and i is a point of this 
manifold. Generally, (ff" is a collection of fields, and then i is a set containing 
also the indices labelling these fields. The hat designates an operator. The 
is an operator in a Hilbert space which is not granted. The workers have 
to build it with their own hands as a representation of the algebra of i^S's. 
For simplicity, (p^ will be assumed boson and real (self-adjoint) but otherwise 
arbitrary. 

The starting point is an operator equation for 

5,((^) + J, = (1.2) 

which is understood as an expansion. It is meant that there is a c-number 
function S'i((p) understood as a collection of its Taylor coefficients at some 
c-number point of configuration space: 

'^»(^) = E -^S^n-i^m^p - cY^ . . . (^ - cy- , (1.3) 

n=0 

and one replaces ip^ in this expansion with an operator. Which c-number field 
d' will be used for this expansion does not matter because it will always sum 
with the operator {ip~cy to make the full quantum field. The expansion point 
c-' is often called "background field" , and there has been much emphasis on it. 
In fact it is completely immaterial. I shall never make this expansion explicitly 
but I shall keep explicit the c-number term of the equation: a source J,;. 
Important are only the following three points. 

(1) The function Si{ip) is local, i.e., it depends only on (p and its finite-order 
derivatives at the point i. 

(2) The function Si{Lp) is a gradient: 



(1.4) 
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i.e., there exists an action S{(p) generating the operator field equations. 
For its derivatives the following notation will be used: 

^--(^)-;^---^^(^)- (^•^) 

Of course, only the total action matters: 

Stot = S{ip) + ifi'Ji . (1.6) 

(3) There is a special condition on the matrix of second derivatives of <?(</?). I 
shall refer to this continuous matrix as S^'- 

Sij{^) = S2{^) . (1.7) 

By locality, 5*2 is the kernel of some differential operator on the base mani- 
fold for which I shall use the same notation 82- It is required that ^2 admit 
a well-posed Cauchy problem in which case it has the unique advanced and 
retarded inverses (Green's functions) G+ and G~: 

SijG^^'' = , = G-''^ . (1.8) 

Because S2 is symmetric, the advanced inverse is the transpose of retarded. 

One may think of S2 as of a second-order hyperbolic operator which it will 

in fact be below but the scheme is more general. It is formalism-insensitive. 
One's field equations may have the second-order differential form or the first- 
order differential form, - the scheme will work anyway. The importance of 
the operator S2 is in the fact that it determines the linear term of the field 
equations and, therefore, governs the iteration procedures. Commute with 
the field equations. Obtained will be a linear homogeneous equation for the 
commutator [ip^\ ip^. Consider the respective inhomogeneous equation and its 
two iterative solutions: one with the advanced inverse for ^2 and the other one 
with retarded. The equation for the commutator is solved by their difference: 

[ip\ ^1 = ih {G+^^{c) - G-^^'(c)) + 0{ip - c) . (1.9) 

In this way the algebra of (p's is built as an operator expansion. This is the 
quantization postulate. 

By the setting of its Cauchy problem, the operator ^2 introduces the con- 
cept of causality. If S2 is a second-order hyperbolic operator, this is the usual 
relativistic causality. But in any case the base manifold will be foliated with 
the Cauchy surfaces of the operator ^2. They will be denoted as U. 

A function of (p that involves (f on only one Cauchy surface 

Q{0) = Q{^\^) (1.10) 

will be called local observable. A state defined as an eigenstate of local ob- 
servables 
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) = ) (1-11) 

will be called local state. This latter name may be confusing because the state 
is, of course, a global concept, and I am using the Heisenberg picture. But the 
local state is associated with a given E: 

\) = \E,q). (1.12) 

Of course, for it to be defined, one needs a complete set of commuting local 
observablcs. I call the Q's obscrvablcs but they may not even be Hcrmitian. 
And I shall consider them linear in (p. If they are nonlinear, I shall make a 
local reparametrization of the field variables so as to make them linear. 

In fact, if one has a complete set of commuting local observables, one has 
already built a Hilbert space. A linear combination 

\E)= j dq^{q)\S,q) (1.13) 

is also a local state associated with E provided that the function !Z'(g) is 
external, i.e., independent of the quantum field t^'. 

Our goal is to learn how to calculate expectation values of field observables 
in a local state, and I shall concentrate on the expectation value 

mm ■ (1-14) 

However, we shall save the effort if we consider another problem first. Namely, 
let us recall what would we do in the case of two local states associated with 
different Cauchy surfaces: 

|ri,gi) = |l), \E2,q2) = \2), (1.15) 

E2> El. 

Here and below, "greater" is a notation for "later". 
The Quantum Boundary- Value Problem 

In the problem where given are two local states (1.15), the field's expectation 
value is replaced with the scalar product 

(Lie) 

which I shall call mean field although it is not mean in any state. 

If our goal was the scalar product (1.16), we would use the Schwinger 
principle 

5(2|1) =i(2|55tot|l) or zero (1.17) 
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whose meaning is this. Consider a variation in the Taylor coefficients of the 
field equations, i.e., in the functional form of the total action. The solution 
for (p^ will respond and will induce a change in the functions Q((p) which will 
induce a change in their eigenstates, and finally there will be a change in the 
amplitude (2|1) induced by a change in the action. The Taylor coefficients are 
local. They can be varied in the region between E\ and Si or outside this 
region. The Schwinger principle (1.17) says that, if they are varied outside, 
the variation of the amplitude is zero. Otherwise, this variation is expressed 
through the variation of the action by (1.17). 

The Schwinger principle is a consequence of the commutation relations but 
it can also be taken for the first principle because one does not need anything 
else. For many purposes (but not all) it suffices to use a specific case of (1.17): 
a freedom of varying the source J. The result of this use is 

_± ^/2|i\ = / (215^ {^'^ ■ ■ ■ ^^") if ^2 > Ji, . . . Jn > ^1 , 

JiJji SiJj \ 0, otherwise . 

(1-18) 

Here T orders the operators (p'', k G Sk, chronologically, i.e., places them in 
the order of following of their E^, and the arrow over T points the direction 
of growth of the time S. 

Let us come back to the operator field equations. Since all (p's in these 
equations are at the same point, one can formally insert in (1.2) the sign of 
chronological ordering: 

"TS.iifi) + J, = . (1.19) 

One may worry about additional terms in (1.19) stemming from the distinction 
between the chronological and ordinary operator products, and the noncom- 
mutativity of ^ with the derivatives in the Taylor coefficients of the equations. 
Because the operators in the products are at the same point, these terms are 
ambiguous expressions whose handling depends on the formalisms and pro- 
cedures used. There is always a happy end: these terms cancel and help to 
cancel similar terms appearing in the subsequent calculations. Therefore, it 
makes sense to use such formalisms and procedures that these terms do not 
appear at all. This is the approach that I shall follow. 

Sandwiching the equation (1.19) between the states (2| and |1), and using 
(1.18), one obtains the following equation for the amplitude: 

(5,(^)+Ji)(2|l) = 0. (1.20) 

Multiply it from the left with (2|1)~^ and pull the factors (2|1) in the argument 
of Si using the fact that this is a unitary transformation: 

[s^(^{2\l)-'^{2\l)^+J?jl = 0. (1.21) 

In the argument, commute the operators: 
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^ + ^)+J.)l = (1.22) 



and use that by (1.18) 

51n(2|l) 



= (/) • (1-23) 



SUk 

The result is the following equation for the mean field: 

Si (^{f) + j^) + ■^^) 1 = • (1-24) 

Equation (1.24) differs from the classical field equation by the operator 
addition 5/5\J to (93). When this operator addition acts on 1, its effect is 
zero, but it will act also on {ip) because the summands ((/j) and 6/5iJ do 
not commute. Where in (1.24) is the Planck constant? It is easy to see by 
dimension that U is just in front of 5/5\J . Therefore, if one wants to expand 
the equations in h, one should expand them in S/SiJ. 

The problem boils down to expanding a function f{A + B) in B when A 
and B do not commute. It suffices to expand the exponential function since 
one can write 

or, equivalently, 

f{A + B)^ e(^+^)'^/'^"/(x) . (1.26) 

x=0 

For the exponential function one has the identity 

e(A+B)x ^ gAx ^.Ay^^(,A+B)y^ ^^ 27) 

which makes the expansion possible. This all works well if the series of com- 
mutators 

e-^Be^ = B+[B,A] + ^[[B,A],A] + ^[[[B,A],A],A]+--- (1.28) 

terminates somewhere as in our case. Indeed, if {(f) = A and 5/5iJ = B, then 

[[B,A],A]=0. (1.29) 
Under condition (1.29) one obtains for an arbitrary function: 

f{A + B)= f{A) + f\A)B + \f"{A)[B, A] + 0{B') . (1.30) 



As compared to the ordinary Taylor expansion, there are several additional 
terms with commutators at each order. 
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A use of the result above in equation (1.24) gives 

Sim) + Isijkm)^-^ + o{h') = -Ji , (1.31) 

Si3m)^-^=-^' + 0{n). (1.32) 

Here the second equation is obtained by differentiating the first one, and it teUs 
us what is 5{ip)/5J. Up to 0{h), it is some Green's function of the operator 
82- Denote this Green's function as 

= + 0{n) . (1.33) 

dJk 

One can work to any order but I shall stop here. We obtain closed equations 
for the mean field: 

+ y^,kM)Gi''m) + 0{h^) = -Ji , (1.34) 

5.,((^))G^'=((^)) = -5f . (1.35) 
The second term in (1.34) is the loop 

Si{{ip))+'-Q+0{n^) = -Ji, (1.36) 

all elements of the loop being functions of {(p). But two questions remain to 
be answered: 

(i) Which Green's hmction is G? 

(ii) What are the boundary conditions to the mean-field equations? 

The answers are again in the Schwinger principle. Equation (1.18) tells us 
what are G and (tp): 

Multiply these expressions by the coefficients that make the linear Q out of 

ip: 

Q{^) = kj0^ , (1.39) 

and send j either to Si or to S2- By the definition of the states |1) and |2), 
one obtains 

Q{{<P) ^ ) = 91 , Qiiv) ^)=Q2, (1-40) 
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jk 



= 



jk 



= . 



(1.41) 



Prom (1.37) it follows also that 



QOk _ Qkj 



(1.42) 



The Green's function G is symmetric and completely determined by the 
boundary conditions (1.41). This completes the determination of the mean- 
field equations (1.34), and for these equations one arrives at a boundary- 
value problem with the boundary conditions (1.40). As a result, the quantum 
boundary-value problem is reduced to a c-number boundary-value problem. 
I say "c-number" rather than "classical" because there are differences, and 
one is the presence of terms 0{h) in the equations, but, as far as the setting 
of the problem is concerned, there is no difference. One arrives at the same 
boundary- value problem for the observable field as in the case of the classical 
states. 

Note that the Green's function G and, thereby, the mean-field equations do 
not depend on the eigenvalues q. The eigenvalues appear only in the boundary 
conditions to the equations. However, G depends on the choice of the observ- 
ables Q themselves and, through them, on the choice of the states |1) and |2). 
Therefore, the mean- field equations are state-dependent. 

Although the Green's function G depends on the choice of the states, 
it possesses two universal properties. One has already been mentioned: G 
is always symmetric. The other one is this. Let us make a variation in the 
operator ^2 and find out how does G respond: 



To answer this question, one can use the Schwinger principle again. The result 
is the following variational law: 



and this law is universal. It is the same for all boundary-value problems. 

The variational law (1.43) is remarkable. It is characteristic of finite- 
dimensional matrices. If a matrix has a unique inverse, then the inverse obeys 
this law. This law is valid, for example, for the inverse of an elliptic operator, 
i.e., for the Euclidean Green's function. It is valid also for the advanced and 
retarded Green's functions: 



S2G = -1 



S26G = -6S2G , 



5G =? 



SG = GSS2G , 



(1.43) 



SG+ = G+SS2G+ , 



SG- = G-6S2G 



(1.44) 



But it is not valid generally, and, in the case of ^2, it is exceptional. 
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The variational law for G has an important implication. Namely, let us 
differentiate the left-hand side of the mean-field equations 

ri{^) = S,{^) + ^5„„„(y.)G""(v5) + 0{n'') (1.45) 



to see if the result is symmetric. One obtains 

= Q + 0{n^). (1.46) 



V V ~2i 



This means that Fi^ip) is a gradient, i.e., there exists an action generating the 
mean- field equations: 

r.(.) = ^. (1.47) 

There is another way to arrive at the same conclusion. Consider a function of 
the mean field defined by the Legendre transformation 

r((v5)) = iln(2|l)-(v5'=)Jfc (1.48) 

where J is to be expressed through (ip) by solving equation (1.23). It is easy 
to see that this function satisfies the equation 

and, therefore, its gradient is the left-hand side of the mean-field equations. 
r{ip) is the effective action. Up to it is of the form 

r{ip) = S{(p) + ^ lndetG((p) + 0{h^) (1.50) 
where the second term is the loop without external lines: 

r{^) = s{^)+(^+o{n^). (1.51) 

The effective action exists for any boundary- value problem but these actions 
are different for different such problems. Only in the classical approximation, 
the action and the equations arc independent of the boundary conditions. 
Let us go over to expectation values. 



The Quantum Initial- Value Problem 



In this problem, given is only one local state (which I shall assume normalized). 
Since the field operators are now sandwiched between the states associated 
with one and the same S: 



10 ViUcovisky 



(1|(...)|1), (1|1) = 1 (1.52) 

one cannot apply the Schwinger principle; there is no room for varying the 
source. One can create this room artificially by inserting a complete set of 
states associated with some later S: 

(l|l)=5^(l|2g)(2g|l), (1.53) 



S2 > Si 

but this alone will not help because the source is varied in both amplitudes, 
and these variations cancel. It will help only if the two amplitudes in (1.53) 
are functions of different sources, i.e., if, instead of (1.53), one introduces a 
function of two independent sources, J and J*: 

z{j*,j) = Y,mj'{H^)j ■ (1-54) 



This amounts to considering two copies of the quantum field: one with the 
source J, the other one with the source J*, and using in (1.54) the amplitudes 
of both. Then one can vary only one source and, after that, make the sources 
coincident. Using the Schwinger principle, one obtains 



(5"Z(J*, J) 



i5iJjj • • • 5i Jj„ 



{1\'T {ifP^ ...ifP^)\l) . (1.55) 



J' = .T 



In this way the expectation values can be calculated. 

The technique of two sources is called time-loop formalism because in 
expression (1.54) one goes forward in time, from Ei to some E2, and then 
back from S2 to Si but with another copy of the quantum field. 

For every partial amplitude in (1.54) we have equation (1.20) 

S^|^-^^+J?j{2q\l)J = 0. (1.56) 

Since the other amplitude in (1.54) does not depend on J, we can linearly 
combine equations (1.56) to obtain 

S^(^-^^+J^^Z{J*,J)=0. (1.57) 

Only one source is active in this differential equation. The other one is a pa- 
rameter. Therefore, we can just repeat the consideration above with Z{J*, J) 
in place of (2|1), and in this way derive the mean-field equations. We obtain 
the loop expansion of exactly the same form as before: 

Sim) + j^Sijkm)G'''m) + oin') = - j, , (i.58) 
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S,,((^))G^'=((^)) = -<5f , (1.59) 

and in these loops we must make the sources coincident. There are only two 
elements in all loops, {^p) and G. Upon setting J* — J, (i^) becomes the 
genuine expectation value 



(/> = 



5\nZ{J*,J) 



SiJk 



and the matrix G is given by the expression 



0{h) 



(52 lnZ(J*, J) 



SiJjSiJk 



(1.60) 



(1.61) 



J* = ,7 



I am using for it the same letter G but it is now a different Green's function 
of the operator £'2. Equations (1.58) with this Green's function in all loops 
are the expectation-value equations. 

The solution of the expectation- value equations is specified completely by 
the initial conditions on Si following from (1.60) but it is not easy to write 
these conditions down in the general terms. Only half of them is obvious: the 
Q's on El are given. To obtain the other half, one would need to find the 
variables canonically conjugate to Q's and calculate their expectation values 
on SiE The same concerns the specification of the Green's function G. This 
issue will be considered in the next lecture where a different approach to it 
will be used. 

Let us consider the state- independent properties of G. First, as seen from 
(1.61), G is symmetric for any initial-value problem: 



(1.62) 



Second, one can apply the Schwinger principle to derive the variational law 
for G. At this point, the initial-value problem differs significantly from the 
boundary-value problem. When the operator 5*2 is varied in the generating 
function (1.54), one can no longer play with only one source because S2 is the 
same for both copies of the quantum field, and, therefore, both amplitudes in 
(1.54) respond. As a consequence, all four matrices of second derivatives are 
generally involved: 



^2 InZ 



5^ \nZ 



6iJj6Uk ' SiJ*6iJl 



6^ InZ 
SUfSLh 



6^ \nZ 
SiJjSiJ^ 



(1.63) 



^ Let Q's be Hermitian, and let P's have c-number commutators with Q's: 
[P, Q] = i. Then the expectation values in the state (1.13) satisfy the initial con- 



ditions 



(Q 



) 



dq^{q)q1r[q) 



) 



— d 

dq^{q) — 1'{q) 



where the overline means complex conjugation. If both Q(<^) and P(<^) are linear, 
these are initial conditions directly for (<^). 
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i.e., the Green's function G^'^, its complex conjugate, and two Wightman func- 
tions: (1|(^^(^'^|1) and its transpose. The Wightman functions can be expressed 
through G^*^ and the advanced or retarded Green's function: 

i(l|(^V''|l) - i(<P')(/) = G^'' - 0+^* + 0(^1) = G''^ - G-''^ + 0{h) . (1.64) 

The result of the calculation is the following variational law for G: 

5G = G-5S2G + GSS2G+ - G-SS2G+ . (1.65) 

It is no more the simple law (1-43) but it is, nevertheless, universal because 
G+ and G~ are state-independent. The variational law (1.65) is valid for any 
initial-value problem. 

The left-hand side of the expectation-value equations has the form (1.45) 
as before but, since the variational law for G is different, the former inference 
about the symmetry of SFi/Sip^ needs to be revised. This inference is no longer 
valid. The advanced and retarded Green's functions arrange it so that 

= when i < j (1.66) 

and 

^^^9^0 wheni>j. (1.67) 

It follows that there is no action generating the expectation-value equations. 

The nonexistence of an action for the initial-value problem is seen also 
from the consideration of the Legendre transform of the generating function 
(1.54). It is now a function of two fields: 

r(v?*, v) = Y In Z{J*, J)-<fiJ + if* J* (1.68) 

where 

6\nZ{J*,J) , ^lnZ(J*,J) 

^= ^L^' ^= 6U^- 

The expectation-value equations are obtained as 



^=(1^1) : 
and, therefore. 

This is not a gradient. 



= -Ji , (1.70) 



(1.71) 
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2 The In- Vacuum State and Schwinger— Keldysh 
diagrams 

Specification of The State 

In order to proceed, I need to specify the state. This wiU be done in several 
steps. 

Step 1. It will be assumed that S2 is a second-order hyperbolic operator, and 
the energy-momentum tensor of the field of small disturbances S^p'^ with the 
action 



S'y(5(/9*V (2-1) 



1 

2 

satisfies the dominant energy condition. 

Step 2. The initial- value surface will be shifted to the remote past: 



. (2.2) 
Consider the operator field equations (1.2)-(1.3): 



00 

Ji+ Si{c)+ Sij{c){^- cy + Y,-,Sin---iM{^- . . .{0- cy- = o . (2.3) 

If c' is some classical solution: 

Si{c) = -Ji, (2.4) 
and (j)^ is an operator solution of ^2 against the background d: 

Si^{cW=Q, (2.5) 

then the field 

= c' + ^\ ieS^-oo (2.6) 

solves the operator dynamical equations asymptotically in the remote past. 
It is a property of ^2 that its solution with smooth data having a compact 
support or decreasing at the spatial infinity decreases also in the timelike 
directions. Then, as i G ^ — oo, the nonlinear terms in (2.3) decrease even 
faster and are negligible. Thus, to build a Hilbert space of states, it suffices 
to build a representation of the algebra of ^'s. 

Step 3. A Fock space will be built associated with the linear field (ji*. This 
amounts to expanding in some basis of solutions of S2{c): 

S2ic)XA = 0, (2.7) 

<^' = x^ai„^ + x^ai^ (2.8) 
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where the overUne means complex conjugation, and the basis functions Xa 
are normaUzed with the aid of the inner product: 

iXA, Xb) = , (xai Xb) = Sab , (2.9) 



(01, (/.2) = -ij (l)iW^(j)2 dS^ . (2.10) 
i: 

Here is the Wronskian of ^2. In this way, the concept is introduccxi of 
some particles detectable in the past. What kind of particles are these, i.e., 
what kind of detectors detect these particles - depends on the choice of the 
basis of solutions but, in any case, the following functions will be chosen for 
the local observablcs Q: 

Q^i^l^) = -^S^"" j XBW^iip - c) dS^ , (2.11) 

s 

S — oo . 

One needs these observables only on the initial- value surface, and, there, they 
coincide with the annihilation operators of the introduced particles: 



Q^iif ) = Oin^ . (2.12) 

The choice of the quantum state will be made in favour of the zero-eigenvalue 
eigenstate of these observables: 

ain^|l) = 0. (2.13) 

This is the vacuum of the introduced particles. 

It follows from (2.6) and (2.8) that the field's expectation value in the state 
(2.13), when taken in the remote past, coincides with the classical solution c*: 

{l\ip^\l)=c\ ieS^-oD. (2.14) 

The ad hoc classical solution c' can then be eliminated completely both from 
the asymptotic form of the quantum field 

(fi' = {ifi') + 4)\ iGE^-Qo (2.15) 

and from the equation defining the Fock modes 

%((V5))(^^' = , ieS^-oo. (2.16) 

Only the mean field itself figures as a background. 

The specification of the state is, however, not completed because the mean 
field in the past remains an arbitrary classical solution: 

Si{{^)) = -Ji, ieS^-oo (2.17) 

and the state itself remains the vacuum of undefined particles. To make the 
final determination, one more step is needed. 



A 
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Step 4- The final choice of the state assumes one more limitation on the 

original action. Namely, it will be assumed that the external source Ji and 
all the external fields that may be present in the action S are asymptotically 
static in the past. This means that, asymptotically in the past, there exists a 
vector field ^'^ such that it is nowhere tangent to any of the Cauchy surfaces, 
and the Lie derivative in the direction of ^'^ of all external fields is zero. 
Specifically, 

C^Ji = , ieS ^ -oo . (2.18) 

If this limitation is fulfilled, then, among the solutions of (2.17) for the 
mean field in the past, there is the static one: 

£j((^»)=0, iGE^-oo. (2.19) 

Choose it. Next, use the fact that, with this choice, the operator S2{{'fi)) 
commutes with the Lie derivative, and choose for the basis solutions of 5*2 ((</?)) 
the functions that, asymptotically in the past, are eigenfunctions of the Lie 
derivative: 

i£^XA = ^ax'a , £.4 > , i G r ^ -00 . (2.20) 

This fixes both the initial conditions for the mean field and the type of particles 
whose vacuum is the chosen state. These are particles with definite energies. 

Since 5*2 is a second-order hyperbolic operator, it contains some tensor 
field, gf^", contracting the second derivatives. The inverse matrix, g^i,, can 
serve and does serve in every respect as a metric on the base manifold. The 
metric enters the original action 5* either as a part of the quantum field ip^ or as 
an external field. In both cases it is subject to equation (2.19). When applied 
to the metric, this is the Killing equation. Thus, we assume the existence, 
asymptotically in the past, of a timelike Killing vector 

The specification of the quantum initial data is now completed. The no- 
tation for the state defined above is 

|l) = |invac), (2.21) 

and its full name is relative standard in- vacuum state. It is "relative" because 
it is relative to the background generated by an asymptotically static source. 
It is "standard" because it refers to the standard concept of particles. It is 
"in" because these particles are incoming. And it is "vacuum" because these 
particles are absent. 

The state should not necessarily be chosen as the zero-eigenvalue eigen- 
state. Since the expectation- value equations do not depend on the eigenvalues, 
they will have the same form for any eigenstate of the annihilation operators, 
i.e., for any coherent state 

ai„^|in a) = a^|in a) . (2.22) 
Only the initial conditions for the mean field will be different: 
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(a iii|(^*|in a) = d + x\(^^ + Xa"''^ > i S ^ -oo . (2.23) 

In addition to the static background d generated by a source, the mean field 
in the past contains now the incoming wave of an arbitrary profile. This is 
the general setting of the classical evolution problem for an observable field 
like the electromagnetic or gravitational field. The fact that the nature of the 
state has changed from classical to quantum did not aff'ect this setting. 

It will be useful to keep comparing the initial-value problem with the 
boundary-value problem. In the latter case, one can define similarly the out- 
vacuum state and specify the quantum boundary data as 

|1) = |in vac) , |2) = |out vac) . (2.24) 
Perturbation Theory 

With this specification of the states, let us come back to the mean-field equa- 
tions. There remains to be obtained the Green's function G(}p) that figures in 
the loops. We need it for an arbitrary background but we have a variational 
law, (1.43) or (1.65), which may be regarded as a difi^erential equation for G'((^) 
with respect to The only thing that is missing and that depends on the 
choice of states is the initial condition to this equation. It suffices, therefore, 
to know G for only one background. 

Then let us do the simplest: perturbation theory around the trivial back- 
ground. A second-order hyperbolic operator with the trivial background is the 
D'Alembert operator with flat metric. Do: 

52((p) = ao + P. (2.25) 

The remainder is a perturbation P. 

In the case of the boundary- value problem, the variational law is (1.43), 
and, therefore, the expansion of G(</?) is of the form 

G{^p) = Go + GoPGo + GoPGoPGq + ... (2.26) 

where Go is G for the trivial background. This expansion is to be inserted in 
the loop in the mean-fleld equations 

^S,MG^\'P)='-Q). (2.27) 

Let for simplicity P be a potential. One obtains the loop expanded in powers 
of P: 

= j dy^... dyn F(x|2/i, . . . 2/„)P(2/i) . . . P(2/„) . (2.28) 

The coefiicients F will be called formfactors. The formfactors are loop dia- 
grams 
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F{x\y) = 



(2.29) 



F{x\yi,y2) 




yi 



2/2 



(2.30) 



with the same propagator for all lines: the trivial-background Green's function 

= Go. (2.31) 

What is Go? With the trivial background and the standard in- and out- vac- 
uum states, it is the Feynman Green's function: 

Go — Gfeynman • (2.32) 

Let us do the same thing for the initial-value problem. The loop in the 

expectation-value equations will, in the same way, be expanded in powers of 
the perturbation, and the expansion will have the same form (2.28), but the 
formfactors will be different because the variational law for G is different. It is 
now (1.65) rather than (1.43). Using this law, one obtains for the formfactors 
three diagrams in place of one: 




(2.33) 



five diagrams in place of one: 



F{x\yi,y2) = 




(2.34) 



and so on. There are two types of propagators in these diagrams: the trivial- 
background G, and the trivial-background retarded or advanced Green's func- 
tion. Respectively, there are two types of lines: 

= Go , = Go or G+ . (2.35) 

In the latter case, the arrow points the direction of growth of time. And what 
is now Go? In terms of the linear field (2.5) it is 
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IgI'' = (in vac\'T{^(j>'')\m vac) (2.36) 

i trivial background 

and differs from the previous case in that the "(out vac|" is replaced by the 
"(in vac|". But, with the trivial background, the vacuum for the linear field 
is stable. The out-vacuum coincides with the in-vacuum. Therefore, 

Go = Gfeynman (again!) . (2.37) 

The diagrams above are called Schwinger-Keldysh diagrams. There is not 
more than one Feynman propagator in every diagram. The remaining ones 
are the retarded and advanced Green's functions organized in a special way 
and with special signs of the diagrams themselves. There is a mystery in 
this special arrangement. What do these diagrams want to tell us? We must 
disclose their secret because working with them directly is not what can be 
recommended. 



Mystery of The Schwinger— Keldysh Diagrams 

One thing is obvious right away. In the diagrams above, there is always a chain 
of retarded Green's functions connecting a given point y with the observation 
point X. Therefore, the formfactor vanishes if at least one of the y's is in the 
future of X. This is the retardation property 

F{x\yi,...yn) =0 when ym > x , Vm . (2.38) 

But this is true of every Schwinger-Kcldysh diagram, and why do they appear 
in the special combinations? What is the role of the Feynman propagator? 

Let us make a Fourier transformation of the formfactor with respect to 
the differences {x — ym) in the Minkowski coordinates: 

/n 
dki...dkn exp(^i ^ km{x - ym)j f{ki, ■ ■ -kn) ■ (2.39) 
m=l 

How come that F possesses the retardation property? It is only that / should 
admit an analytic continuation to the upper half-plane in the timelikc com- 
ponents of k's. Then, for y^ later than x, we shall be able to close the inte- 
gration contour in the upper half-plane of k^, and the integral will vanish. 
There should be a function of complex momenta /(zi, . . . Zn) analytic in the 
upper half-planes of and such that f{ki, . . . kn) is its limiting value on the 
real axes: 

/(fci,...fc„) = /('2l,---^n) _ ,0 , . • (2-40) 

Let us build this function. 

All diagrams in a given-order formfactor are similar. They all arc integrals 
over the momentum circulating in the loop, and the integrands are identical. 
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The difference is only in the integration contours. Thus any diagram in the 
lowest-order formfactor f{k) is of the form 

/■ polynomial in momenta 

c 

There are, generally, as many factors in the denominator as there are propa- 
gators in the loop, and each factor contains two poles. The contour C passes 
round them in accordance with the type of the propagator. One of the three 
rules applies to each pair of poles: 

retardation rule, 

advancement rule, 

Feynman rule. 



Let us now shift the external momentum k'^ to the complex plane. The 
poles will shift to the complex plane but we shall also deform smoothly the 
contour so that it do not cross the poles. In this way one can build a function 
of complex momenta for each Schwinger-Keldysh diagram. Thus the lowest- 
order formfactor with complex momentum, f{z), is a sum of three functions: 

f{z) = JdpJ dp°i...)+ jdpj dp\...)- jdpj dp\...), (2.42) 
Ci C2 C3 

and the contours Ci, C2, C3 for in the upper half-plane are shown in Fig. 1. 
By considering the pinch conditions, i.e., the conditions that the poles pinch 
the integration contour, one can check in each case that these functions can 
have singularities only on the real axis. Therefore, if we consider them in the 
upper half-plane, they are analytic, and their limits on the real axis are our 
original diagrams. 

There remains to be understood what are these functions. Since the inte- 
grands are identical, the sum of the integrals in (2.42) is the integral over the 
sum of the contours 

f{z) = j dp J dp° (...). (2.43) 
C1+C2-C3 

Sum up the three contours in Fig. 1. The resultant contour is such that every 
pair of poles is passed round by the Feynman rule. It may be called Feynman 
contour. 

But the Feynman contour defines also the in-out formfactor (2.29) in which 
both propagators are Feynman, except that the in-out formfactor is not the 



20 Vilkovisky 



+ 



+ 



p° plane 



TV. 



C2 
— 



Tv. 



SUM: 




--FEYNMAN 



Fig. 1. Integration contours for the three diagrams in the lowest-order formfac- 
tor (2.42). The sum of the contours is the Feynman contour. 



limit of 7(2;) from the upper half-plane. It is this limit on only half of the 
real axis, and on the other half it is the limit from the lower half-plane. The 
in-in and in-out formfactors are different boundary values of the same complex 
function having a cut on the real axis: 

in-in: /W = ^ ^0 + ' (2-44) 

in-out : /W = /WLo = (i + i,)fcO' (2-45) 
and the function itself is the integral over the Feynman contour 

/(z) = y dp y (...). (2.46) 

CpEYKMAN 

The same is true of all n-th order formfactors, and this is a disclosure of 
the mystery. In each case, the set of Schwinger-Keldysh diagrams is just a 
splitting of one Feynman diagram whose purpose is to display the retardation 
property and in this way to tell us which boundary value is to be taken. 
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Reduction to The Euclidean Effective Action 

The Feynman contour is famous for the fact that, when the external momenta 
arc on the imaginary axis, the Feynman contour is the imaginary axis itself. 
With all the momenta imaginary, both the external ones and the one circulat- 
ing in the loop, this is the Euclidean formfactor. Then we can start with the 
calculation of the Euclidean formfactor and next analytically continue it in 
momenta from the imaginary axis to the real axis either in the way shown in 
Fig. 2(a) or in the way shown in Fig. 2(b). In the first case we shall obtain the 
in-out formfactor, and in the second case the in-in formfactor of Lorentzian 
theory. It is invaluable that loops can be calculated Euclidean. 




Then let us make one more step. A formfactor with the Euclidean momen- 
tum can be put in the spectral form 

/(fc) = / dm^ ^'f^ + a polynomial in fc^ , (2.47) 



fc2>0 

with some spectral weight p{m^), the resolvent l/(m^ + fc^), and a polynomial 
accounting for a possible growth of f{k) at A;^ ^ oo. There are similar forms 
for the higher-order formfactors. If the formfactor is in the spectral form, 
the procedure of analytic continuation boils down merely to replacing the 
Euclidean resolvent with the retarded or Feynman resolvent: 



oc 

dm^ — 7, 777^ — 72 + ^ polynomial in 





W? - {k° + i£)2 + 

(2.48) 
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. ,, . f , 9 aim?) , . , . , 2 

m-out : j[k) = am — — ^ ; — h a polynomial m k . 





— k^"^ + k — ie 

(2.49) 

Note that the spectral weight is the same in all cases: the one of the Euclidean 
loop. Thus, the problem boils down to obtaining the spectral weights of the 
Euclidean formfactors. 

Then back from the Fourier-transformed formfactors to the formfactors 
themselves, and from the formfactors to the mean-field equations. For the 
loop in these equations expanded in powers of the perturbation, we obtain an 
expression of the following form: 



{ci+C2ao + ...)P{x) 

DO 

1 



+ / dm' pim')^^^^^—^P{x) 







+ 



J dm\dm\dm\p{m\,m\,m^l) 

n4] ( -T^nx) 



mi - Dn 



l2 - ^0 



□n ^ ' \mi-U 



+ ... . (2.50) 

Here the first term is local. It comes from the polynomial in the spectral 
form. The remaining terms are nonlocal but expressed through the resolvent 
which is a Green's function of the massive operator Dq — m^. It is initially 
the Euclidean Green's function since we are calculating the Euclidean loop. 
For the Lorentzian equations, we arrive at the following rule. To obtain the 
expectation-value equations in the in- vacuum state, replace all the Euclidean 
resolvents in (2.50) with the retarded Green's functions. To obtain the mean- 
field equations for the in-out problem, replace all the Euclidean resolvents 
with the Feynman Green's functions: 

< Euclidean, 
Retarded, (2-51) 
Feynman. 

At every level of expectation-value theory, there are proofs that the 
expectation-value equations possess two basic properties: they are real and 
causal. Causality is the retardation property discussed above. But it is not 
enough to have proofs. These properties should be manifestly built into the 
working formalism. Expression (2.50) offers such a formalism. Since the re- 
tarded resolvent secures the causality and is real, this expression is manifestly 
real and causal. 
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But even this is not enough. The theory may possess symmetries, and 
one may want these symmetries to be manifest. To this end it will be noted 
that, although expansion (2.50) is obtained in terms of the trivial-background 
resolvent l/{m? — Do), it can be regrouped so as to restore the full-background 
resolvent ^ ^ 

m2-52 " m2 - Do - P ^^-^^^ 
at each order. It docs not matter whether this regrouping will be made in 
the expectation-value equations or in the Euclidean equations because the re- 
tarded and Euclidean Green's functions obey the same variational law (1.43): 

1 - 1 . + .... (2.53) 



m? — Do — Si w? — S2 - 

This proves that the rule of replacing resolvents applies to the full-background 
resolvents as well as to the trivial-background ones. The latter fact is im- 
portant because the Euclidean loops can be calculated covariantly from the 
outset, and the transition to the expectation- value equations by replacing 
the full-background resolvents does not break the manifest symmetries. The 
expectation-value equations are obtained in as good an approximation as the 
Euclidean equations are. 

There remains to be made a final observation. For the Euclidean equations, 
there is an effective action: 

because the; variational law for the Euclidean Green's function is (1.43). It is 
invaluable that loops can be calculated without external lines. This reduces 
the calculations greatly, helps to control symmetries, helps to control renor- 
malizations. 

Thus, at the end of the day, we conclude that there is an action that 
generates the expectation- value equations but it does so indirectly, i.e., not 
through the least-action principle. To make this clear, consider (for the illus- 
trative purposes only) any quadratic action: 



r(</.) = i J dxipfiao)ip. 



Whatever the operator /(Dq) is, in the variational derivative it gets sym- 
metrized: 

^ = 1 (/(Do) + r (Do)) <p = f^^^iOo)^ . 

Assuming that the function /(Oo) is in the spectral form 

00 

/(□o) = J dm^p{m\ ^ 



2 
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one obtains the variational equations with the symmetrized resolvent: 



°° sym 



/ 





1 



These cannot be the expectation-value equations since they are not causal. 
But, through the derivation above, we know how to correct this: just to re- 
place the symmetrized resolvent with the retarded resolvent. The corrected 
equations 

ret 



/ 



1 

- Do 



dm p{m^) { 32 — — ] ^ = -J ■ 





do not already follow from any action although indirectly they do. Only if the 
action r{ip) is local, i.e., the function /(Do) is polynomial, the least-action 
principle holds directly. 

Two precepts should be kept in mind when using the formalism above. 
First, the replacement rule concerns the resolvents of the formfactors and not 
the propagators in the loop. The loop should be calculated Euclidean. Hence 

First Precept: first do the loop, next replace the resolvents. 

Second, the replacement of resolvents is to be made in the equations and not 
in the action. It does not make sense to make it in the action. Hence 

Second Precept: first vary the action, next replace the resolvents. 

We thus go over to the calculation of the Euclidean efi^ective action. 



3 The Effective Action 
The Operator S2 

The 1^* is a set of fields for which a more explicit notation will now be used: 

= ^%x) . (3.1) 

The operator S2 acts on a small disturbance of 99* and is a second-order 
differential operator 

SijSif^ = {X^^d^d, + Y^^d^ + Zab) <5/(x) . (3.2) 

The generality of this operator will, however, be restricted by the condition 
that the coefficient of the senior term factorizes as 

Kb='^ab9^'', deta;„6^0, det^'^'^^O. (3.3) 
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In this case, the operator (3.2) is said to be diagonal, or minimal, or nonexotic. 
Condition (3.3) is too restrictive and not necessary. It can be replaced by a 

more general condition 

det{X^^n^n,) = C{g'"'n^n,f Vn^ , d = dima, C 5^ , det^^^'^O, 

(3.4) 

and even this condition can be generalized. Higher-order and first-order oper- 
ators can also be considered but, in all of these cases, the Green's functions 
of 5*2 are expressed through the Green's functions of a diagonal second-order 
operator. The case (3.3) is basic. 

In the case (3.3), the matrix coab can be factored out: 

SijSip^ = ujacm5<^\x) , (3.5) 

and a covariant derivative can be introduced: 

V^^^« = (<5,«5^ + ^^?)V (3.6) 

so as to absorb the first-order term: 

= S^g^-W^W, + PS . (3.7) 

This is the final form of ^2. A short notation will be used: 

H=Di + P (3.8) 

where 

□ = s^-^V^V, , (3.9) 
and the hat designates a matrix in a,b: 

i = S^, P = P,\ tr P = P:, etc. (3.10) 

The matrix Wab may be regarded as a local metric in the space of fields. The 
symmetry of ^2 implies that this matrix is symmetric, covariantly constant, 
and converts P into a symmetric form: 

<^ab = l^ba , V^Wab = , (3.11) 
P>e6 - Pb^ca = . (3.12) 

The dominant energy condition implies that u)ab is positive definite. The ma- 
trix g^"^ is the inverse of the metric on the base manifiold. Since we are con- 
sidering Euclidean theory, this metric is positive definite too. 

Apart from the algebraic factor Wac in (3.5), the operator 5*2 contains three 
background fields: 

5"% V^, P (3.13) 

i.e., the metric, the connec;tion (or covariant derivative), and the matrix poten- 
tial. And where is the original background (p of 5*2 (</?)? When ^2 is calculated 



26 Vilkovisky 



from the action S, the metric, connection, and potential are obtained as func- 
tions of the original set of fields ip, but from now on it does not matter. The 
effective action is expressed in a universal manner through the fields (3.13) 
only. 

The strengths of the fields (3.13) are respectively the Riemann tensor, 
the commutator of covariant derivatives, and the potential which is its own 
strength: 

i?a/3,i. , [V^, v.] =^1,^, p. (3.14) 

I shall call these field strengths curvatures and use for them the collective 
notation 

[Ra|3^^u,n^., P)=^. (3.15) 
The following contractions of the curvatures will be called currents: 

= vn^. , (3.16) 
J^iv = — -j^QiivR , J = g'^'^Jjiv ■ (3-17) 

The currents are conserved: 

V = , V J^. = . (3.18) 

If all the curvatures vanish, the background is trivial. The effective action is 
a functional of the curvatures (3.15). 



Redundancy of The Curvatures 

The effective action is a nonlocal functional of the curvatures, and this fact 
conditions a certain simplification. 

Since the commutator curvature is a commutator, it satisfies the Jacobi 
identity, and so does the Riemann curvature: 

v^Te^. + v^n^^ + = o , (3.19) 

VjRa0nv + ^vRafijn + ^ iJ,Ral3vy = . (3.20) 

Act on these identities with V'. In the first term, the operator □ forms, and 
in the remaining terms commute the covariant derivatives. The commutator 
brings an extra power of the curvature. The equations obtained 

aiZ^, + 0(3?2) = 2V[,J^] , (3.21) 

ORa/3^,. + 0{^^) = 4V[^V(a (| J,]/3> - iff,]^) jj (3.22) 

hold identically and have the form of inhomogeneous wave equations, the role 

of inhomogeneity being played by the currents. In (3.21), (3.22), the brackets 
of both types [ ] and ( ) denote the antisymmetrization in the respective indices. 
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The equations (3.21) and (3.22) are nonlinear but they can be solved by 
iteration. The result is that the commutator and Ricmann curvatures get ex- 
pressed in a nonlocal fashion through their currents and an arbitrary solution 
of the homogeneous wave equation 

□^wave^O^ □i?X' = 0. (3.23) 

If the metric is Lorentzian, this solution is fixed by initial data which can be 
given in the remote past. It follows that the commutator and Riemann cur- 
vatures are specified by giving an incoming wave and the current J. This fact 
underlies the Maxwell and Einstein equations. They fix the currents J. Adding 
initial conditions to these equations specifies the connection and metric. 

In the present case, since the metric is Euclidean, there are no wave solu- 
tions: 

nj;r=o, d;°=o, (3.24) 

and the Green's function !/□ is unique. Therefore, the commutator and Rie- 
mann curvatures are expressed entirely through their currents: 

= ^2V[,J^] + 0(J2) , (3.25) 

Ra0^.u = 54V[^V(„ (^J,,]^) - ^g,.]p)J^ + 0{J^) . (3.26) 

Thus, the curvatures are redundant because there are no waves in Eu- 
clidean theory. Owing to this fact, the set of field strengths (3.15) reduces 
to 

[j^., Jf., P) , (3.27) 
and the effective action is a functional of the reduced set. 

The Axiomatic Effective Action 

To what class of functional does the effective action belong? One can say in 
advance that this should be a functional analytic in the curvature. Indeed, 
the first variational derivative of the effective action taken at the trivial back- 
ground should vanish because, in the absence of an external source, the rel- 
ative vacuum becomes the absolute vacuum. The trivial background should 
solve the mean-field equations in the absolute vacuum. Higher-order varia- 
tional derivatives taken at the trivial background determine the correlation 
functions in the absolute vacuum. They may not vanish but neither should 
they blow up. 

The analjfticity suggests that the effective action can be built as a sum of 
nonlocal invariants of A''-th order in the curvature: 



JV 



(3.28) 



28 Vilkovisky 



Nonlocal invariant is, however, an uncertain concept. Even local invariant of 
A^-th order in the curvature is a concept that needs to be refined but this 
is easy to do. The most general local monomial that can be built out of the 
available quantities yields an invariant of the form 

J dx g^/' (Vi...Vi)(V2...V2)... SRijR2 ...^N + 0[3fi^+i] . (3.29) 



This monomial is a product of N curvatures and k covariant derivatives, all 
indices being contracted by the metric. In (3.29), the labels 1,2,... point out 
which derivative acts on which curvature but all the curvatures are at the same 
point, and the total number of derivatives is finite. Of course, the curvature 
sits also in the covariant derivatives and in the metric that contracts the 
indices. Therefore, the A^-th order invariant can only be defined up to terms 
0[3fJ'^"'"^]. In particular, the covariant derivatives in (3.29) can be commuted 
freely because the contribution of a commutator is already 0[5R^+^]. 

One may now consider a class of nonlocal invariants that can formally be 
represented as infinite series of local invariants: 

/CXD 
dxg^/^ ^Cfc (Vi...Vi)( V2...V2) . . . jRig?2 . . .g?jv + 0[jR^+^] . (3.30) 
fc=o r ' 

k 

Here Ck are some dimensional constants. It can be seen that this is the needed 
clas^. The number of curvatures in (3.30) is N but the number of deriva- 
tives is unlimited. Only a finite number of derivatives can contract with the 
curvatures. The remaining ones can only contract among themselves. If two 
derivatives acting on the same curvature contract, they make a □ operator 
acting on this curvature: 

Vi^^Di, V2' = n2, ... . (3.31) 

If two derivatives acting on different curvatures contract, the contraction can 
again be written in terms of the □ operators: 

2V1V2 - (Vi + V2)' - Vi^ - V2' 

= 01+2 - Di - 02 (3.32) 

but there appears a □ operator acting on the product of two curvatures: 

□ i+2 5Ri5R23?3... = □(5R5R) SRs... ■ (3.33) 

As a result, (3.30) takes the form 



^ To see it, consider any diagram with massive propagators and expand it formally 
in the inverse mass. The method that accomplishes this expansion is known as 
the Schwinger-DeWitt technique. 
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r, 



N 



X (v...3?iV...3?2...V...3?jv) +0[3?^+i] . 

^ V " 

contraction 

(3.34) 

There remains an infinite scries in the □ variables, and these variables them- 
selves are operators acting on the curvatures in a given contraction. The re- 
maining series is some function of the □ variables: 

Fn = f dxg^/^F {0^,02,0^+2, ■■■) (v...3?iV...3?2...V...3?iv) +0[?Si^+^] . 

•' V ' 

contraction 

(3.35) 

This is the general form of a nonlocal invariant of A''-th order in the curvature. 
The function F is a formfactor. 

There is, in addition, the identity 

Vi + V2 + ... + VAr = (3.36) 

which reduces the number of variables in the function F. The sum in (3.36) 
is a derivative acting on the product of all curvatures, i.e., a total derivative. 
Total derivatives vanish because the curvatures may be considered having 
compact supports. Thus invariants of first order in the curvature can only be 
local because any derivative is a total derivative. Therefore, the first-order 
formfactors are constants: 

N=l : F = const. (3.37) 

At the second order, all formfactors are functions of only one argument be- 
cause the remaining arguments can be eliminated by integration by parts: 

N = 2 : F = F{0^) , (3.38) 

□2 = ni, 01+2 = 0. 

At the third order, all formfactors are functions of three individual D's because 
the n's acting on pairs can be eliminated: 

A^ = 3 : F = F(ni,n2,n3) , (3.39) 

□ l+2 = Ds , Ol+S = ^2 , n2+3 = Dl • 

The n's acting on pairs appear beginning with the fourth order in the curva- 
ture and are parameters of the on-shell scattering amplitudes. 
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Nonlocal invariants of a given order make a linear space in which all pos- 
sible contractions of N curvatures and their derivatives make a basis, and 
the formfactors play the role of coefficients of the linear combining. The basis 
can be built by listing all independent contractions. The effective action is an 
expansion in this basis with certain coefficients-formfactors: 



r = Ti + Tii + Tin + . . . 



C1R + C2 trP 



(3.40) 
(3.41) 



Tii = J dxg^/Hr [r^, Fi{n) R^'- 



+ RF2{a) R 

+ PF3{a) R 

+ PF4{D) P 



(3.42) 



Aii = 



tr [ Fi(ni, Da, □3)^1/^2^3 

+ F2(ni, 02, Ds) ni\n2"i3Ti3^^. 
+ ... 

+ F29(ni, 02, 03) VAVCTi?^''V«V/3i?rV;,V.i?^'^ 

(3.43) 

In the first-order action (3.41), there are 2 basis contractions: the Ricci scalar 
and the trace of the matrix potential, and the formfactors are constants. In 
the second-order action, there arc 5 independent contractions listed in (3.42). 
In the third-order action, there are 29 basis contractions, examples of which 
are given in (3.43). Here I shall stop because, for the problems of interest, 
the third order is sufficient. The reason for that will be explained in the next 
lecture. 

In the expressions above, the basis invariants are written in terms of the 
curvatures but they can be rewritten in terms of the conserved currents. Note 
also that the operator arguments of the third-order formfactors F commute 
because they act on different objects. Since the arguments commute, the func- 
tions F themselves are ordinary functions of three variables. 

Thus, even before any calculation, we have an ansatz for the effective 
action, with unknown formfactors. We need them in the spectral forms 



00 

-I 



dm 



2 Pk(m2) 



□ 



-|- a polynomial in □ 



(3.44) 
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OO 

F,(ni, □„ 03) = / drnldmldml , '^^'2"}:"^ , ^ , , (3-45) 
J (mf - □i)(m^ - □2)(to3 - na) 



and then we can proceed directly to the expectation-value equations. Un- 
known are only the spectral weights. These are to be calculated from the loop 
diagrams but there is an alternative approach. One can look for the general 
limitations on the spectral weights stemming from axiomatic theory. These 
limitations may be sufficient to solve one's expectation-value problem. In this 
case, the solution will prove to be independent of the details of the quantum- 
field model and the approximations made in it. Moreover, the effective action 
above does not refer even to quantum field theory. It is an action for the 
observable field, and its implications may be valid irrespective of the under- 
lying fundamental theory. Only certain axiomatic properties of the spectral 
weights may be important. There is an example in which this approach has 
been implemented [53]. 

Here, the axiomatic approach will not be considered. Let us see how the 
effective action is calculated from loops. 



Heat Kernel 

Consider any diagram in the effective action 



(3.46) 



and, for every propagator, write 

00 

=-1= J dse^" . (3.47) 



The kernel of the exponential operator 

e'"S{x,y) = K{x,y\s) (3.48) 

(and the operator itself) is called heat kernel, and the parameter s is often 
called proper time. Both names are matters of history, and a matter of physics 
is the fact that H is negative definite. The matrix P in (3.8) may spoil the 
negativity but, since it is treated perturbatively, as one of the curvatures, this 
does not matter. 

Upon the insertion of (3.47), the diagram remains the same as before but 
with the heat kernels in place of the propagators, and the integrations over 
the proper times will be left for the last: 
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(3.49) 



The one-loop effective action is the functional trace of the heat kernel, inte- 
grated over s: 

(^^= ^Indct-^ = — J dxtvK{x,x\s) . (3.50) 

Thus, one is left with diagrams with the heat kernels. It will be seen in a 
moment why this is better. 

The expansion rule for the exponential operator has already been consid- 
ered in (1.27). There remains to be presented the lowest-order approximation 
for the heat kernel: 

K{x, y\s) = j^-^ (e-<^(-'^)/2^a(a;, y) + 0[K]) , (3.51) 

D = dimension of the base manifold. (3.52) 

At the lowest order in the curvature, the potential P does not affect this 
expression but the metric and connection do. As mentioned above, covariant 
expansions cannot be rigid. In (3.51): 

2a{x, y) = (geodetic distance between x and y)^ (3.53) 

in the metric entering the operator H. The connection entering the operator H 
defines a parallel transport along a line. Parallel transport is a linear mapping, 
so there exists a propagator of parallel transport (the matrix that accomplishes 
this mapping). In (3.51): 

a{x, y) = propagator of the parallel transport from y to a; (3.54) 
along the geodesic connecting y and x. 

The geodesic comes from the metric, and the parallel transport from the 
connection. 

The two-point functions (3.53) and (3.54) are the main elements of the 
Schwinger-DeWitt technique mentioned above and the basic building blocks 
for all Green's functions: of the hyperbolic operator H, and of the elliptic 
operator H, and the heat kernel. What is special about the heat kernel? 
Special is the fact that, as seen from expression (3.51), the heat kernel is 
finite at the coincident points. Green's functions of the hyperbolic and elliptic 
operators are singular, and this is normal. Abnormal is the fact that in the 
loop diagrams they appear at the coincident points. Finiteness of the heat 
kernel at the coincident points is a bonus owing to which all diagrams with 
the heat kernels are finite. 
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The divergences of the loop diagrams reappear in the proper-time integrals 
in (3.49). These integrals diverge at the lower limits. At this stage, one more 
advantage of the heat kernel comes into effect. Namely, the manifold dimension 
D enters only the overall factor in (3.51). Apart from this factor, the expansion 
of the heat kernel in the curvature docs not contain D explicitly. Therefore, 
loops with the heat kernels are calculated once for all dimensions, and then 
the knowledge of the analytic dependence on D enables one to apply the 
dimensional rcgularization to the proper-time integrals. One integrates by 
parts in s keeping Re£) < 4 and next goes over to the limit D ^ A. For 
example, 

oc oo 

/ = 2^^(0) - / + ^ (2 - D/2) . (3.55) 



The dimensional rcgularization annihilates all power divergences. Only the 
logarithmic divergences survive and take the form of poles in dimension. These 
poles affect only the polynomial terms in the spectral representations of the 
formfactors. They appear in the coefficients of the polynomials, thereby mak- 
ing these coefficients indefinite. As a consequence, the local terms of the ef- 
fective action will have indefinite coefficients. I shall come back to this issue. 

After the substitution of the heat kernels for the propagators, the calcu- 
lation of loops becomes an entertaining geometrical exercise. 



Loops and Geometry 

The heat kernel involves a and a. The derivative of a 

X 

Wa{x,y) = (T^'{x,y) ^ ^a^{x,y) (3.56) 

is the vector tangent to the geodesic connecting y and x, directed outwards, 
and normalized to the geodetic distance between y and x: 

g^^a^'a'' = 2a , cr^ =0, detVci^ 7^0. (3.57) 



x=y 



The normalization condition is a closed equation for cr which together with 
the conditions at the coincident points can serve as the definition of a. The 
defining equation for a together with the condition at the coincident points is 

a''V^,a{x,y) =0 , a =1. (3.58) 

The determinant 

det(v;^VXx,y)) =5V2(3,y/2(y)^(3.^y) (3.59) 
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is known as the Van Vleck-Morette determinant. It is responsible, in partic- 
ular, for a caustic of the geodesies emanating from x or y. 

The vector can be used to expand any function in a covariant Taylor 
series. For a scalar, this series is of the form 

f^y^ = E Hi! • • • '^''"^^1 • • • ^^n/w • (3-60) 

ra=0 

If / is not a scalar, it should at first be parallel transported from y to x: 

f{y) = diy, x) ^^"^ ...a^-V^,... V ,J{x) . (3.61) 

n=0 

The covariant Taylor expansion is a regrouping of the ordinary Taylor ex- 
pansion. Whatever the connection is, it cancels in this series. The series can 
formally be written in the exponential form 

f{y) = aiy, x) exp (~a^V^) f{x) (3.62) 

which will be of use below. Two-point functions expanded in this way get 
expressed through their covariant derivatives at the coincident points. Thus 

A{x, y) = l + ^R^.a^'a" + ... . (3.63) 
A loop always involves the ring of a's 



a{x,xi)a{xi,X2) ■ ■ .a{xn,x) , 




(3.64) 



i.e., the parallel transport around a geodetic polygon. The ring of two a's is 
the parallel transport there and back along the same path. Therefore, 

a{x,xi)a{xi,x) = 1 . (3.65) 

The ring of three a's is the parallel transport around the geodetic triangle. It 
involves the commutator curvature, and the curvature terms can be calculated: 

a{x,xi)a{xi,X2)a{x2,x) = i + ^'R.af3<7i°'cr2'^ + ■ ■ ■ , (3.66) 




(3.67) 



This is sufficient because any polygon can be broken into triangles: 
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(3.68) 

Solution of the geodetic triangle is also involved. In the notation of (3.67), 
(a^{xuX2)y =ai^+a2''-2aia2-^R^au0C7i^(Ti''a2"a2'^ + ... . (3.69) 

Here the first two terms make the Pythagorean theorem, the third term ac- 
counts for the angle not being the right angle, and the terms with the Riemann 
curvature can be calculated. 

The above is to give a flavour of what loops imply. 




Calculation of Loops 

The heat kernel calculates loops with a remarkable elegance. As an example, 
consider the contribution of the second order in the curvature to the effective 
action. The respective one-loop diagram contains two curvatures 9? and two 
heat kernels with the proper times si and S2- 

si 

.31 



+ 0[^' 

= J dxg'/^ J dyg'^^^ix)K{x,y\s,)k{x,y\s2Wiy) + 0[^^] . (3.70) 

Suppose that the calculation only needs to be done with accuracy 0[3?^]. Then 
one can insert in (3.70) the lowest-order approximation for the heat kernels. In 
this approximation, the rings of a's collapse to 1, and the remaining a's always 
transport the 3?'s to the same point arranging their complete contraction. 
With the o's and the numerical coeflacients omitted, the diagram (3.70) is of 
the form 



x»(x)exp(-^)exp(-^)9!(„). (3.71) 

But the exponents here simply add, and the two heat kernels turn into one 
with a complicated proper-time argument: 

dxg^'^ j dy5i/2g^(^)exp(^-^-i^a(x,2/)) ^{y) 



(siS2)^/2 
~ (si + 



J dxg'/^J dyg^'^^{x)K[x,y 



S1S2 



Sl + S2 



3?(2/) . (3.72) 
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One only needs to rewrite this heat kernel in the operator form: 

, J ,n/2 I dxg^'^^e^^ ^^v) ' (3-73) 

and the loop is done. The proper-time integral 

oo oo 

5:t^»p (^°) = no) (3.T4) 



is the formfactor. 

What has happened? The propagators in the loop glued together, and the 
loop turned into a tree: 

•<o* — • — • ^ p-^^' 

This is what means to do the loop. It means to turn it into a tree. The role 
of the propagator in the tree is played by the formfactor F{U). 

Consider now any multi-loop diagram with parallel propagators. It turns 
into a tree 




(3.76) 



in a completely similar way. The inverse proper times add: 



1 r . • . — 

Si S2 Stotal 

(the law of parallel conductors). There is nothing to do. 

For more than two curvatures a more powerful method is used. Consider 
the diagram 




(3.77) 



and suppose again that it is needed only up to the next order in the curvature. 
Then, with the a's and the numerical cocfhcionts omitted, it is of the form 



1/2 



X exp ( — ) 3R(x)3?(yi)3?(y2) • (3.78) 
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Choose one of the vertices, say x, to be the observation point of the effective 
Lagrangian. One of the curvatures, is already there. Shift the remaining 

curvatures to x using the covariant Taylor series: 

^{yi) = exp ^{x) , (3.79) 

ai^ =a^'{x,yi) , i = l,2. (3.80) 

Next, consider the geodetic triangle with the same vertices as in the diagram. 
For the geodesies connecting x with y^, write 

2a{x,yi) = {aif , (3.81) 

and, for the geodesic between the y's, use the Pythagorean theorem: 

2a(yi, ya) = (ai)' + (aa)^ - 2aia2 + 0[3?] . (3.82) 

Finally, replace the integration variables: 

yi^^ai^, y2^^cj2^. (3.83) 

The Jacobian 



dyi'' 



i^A-~\x,y,) = ^^{l + Om 



(3.84) 



removes the measure g^/'^ from the integral in yi and brings an extra g^f^ to 
the integral in x. Expression (3.78) takes the form 



(S1S2S3 



L^/dx/'^(,.'^W)7<i,,<i,.exp(-£!-g! 

(71 ''V^^ - (72" V^^) ^{x)^i{x)^2{x) . (3.85) 



CTl^ + fT2^ - 2(Ji(72 

Here the labels 1,2 on V ^ and 3? point out which acts on which 3?. The 
operators W ^ figure as parameters in the integral, and, up to the next order 

in 3?, they commute. Since the parameters commute, the integral in (Ti", 

is an ordinary Gaussian integral. Do it. The extra factor (5^/^(0;)) cancels, 

and the result is 

-B(si,S2,S3) j dxg^^^exp ^ ^ 6ifc(si, S2, S3)ViVfc^ ^{x)^i{x)^2{x) 

(3.86) 

where B{si,S2, S3) is some function of the proper times, and the exponent is 
a quadratic form in Vi, V2 with s-dependent coefficients. The loop is done. 
The integral 
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y dsids2ds3B{si,S2,S3)exp i ^ &ife(si, S2, S3)ViVfe 



i,k=l 

2 V7 2 



= F(Vl^V2^VlV2) (3.87) 

is the formfactor. Integration by parts in x brings it to the □ arguments: 

F{Vi^ V2^ Vi V2) ^ F(Vi2, V2^ V^) . (3.88) 

The effect of the calculation above is again that the loop is turned into a 
tree: 



(3.89) 



The vertex of the tree is the formfactor F(Vi^, V2^, Va^). This method applies 
to any diagram with the heat kernels. One only needs to do Gaussian integrals, 
and the result is always the exponential of a quadratic combination of Vs. 
The formfactor is a function of the products VjVfc. 

The One-Loop Formfactors 

The result of the proper-time integrations depends essentially on the dimen- 
sion D. For D = A, the one-loop formfactors in the effective action (3.40) are 
as follows. 

With one exception, all second-order formfactors are logs: 



Since 



^4(n) = l^^H-^) + const. , (3.93) 



- In(-n) = / dm'^ h const. , (3.95) 

' — □ 
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these expressions have the spectral forms (3.44) with definite spectral weights 
and indefinite additive constants (polynomials of the zeroth power). Respec- 
tively, the effective action contains a set of local terms with unspecified coef- 
ficients: 

r = 7T77-T^ / dx [^i? + ca tr P + c^R^^R'''' + aR^ 

2{ATTy J \ 

-Hc5tr(PF) +C6tr(7^^^7^^'') + -^i?trF+ nonlocal terms) . (3.96) 

18 / 

The nonlocal terms are specified completely. 

The third-order formfactors have no polynomial terms and indefinite coef- 
ficients. The simplest third-order formfactor is i^i(ni, O2, Ds) in (3.43). It has 
the spectral form (3.45), and its spectral weight pi{rn\,'m\,m\) is obtained 
as follows. Consider a triangle of three spectral masses 

A = area of the triangle. 

7713 

It can be built only if every mass is smaller than the sum of the two others. 
The spectral weight pi is zero if the triangle cannot be built. Otherwise, it is 
proportional to the inverse area of this triangle: 

X 9{mi + 7712 — 'mz)0{'mi -f 777,3 ~ m2)d{m2 + — mi) . (3.97) 

The remaining 28 third-order formfactors are expressed through Fi and 
are tabulated [55j . The tables contain various integral representations of the 
formfactors, and their asymptotics. 

The loop of the minimal second-order operator with arbitrary metric, con- 
nection, and potential is called standard loop because every calculation with 
it is done once, and the results can be tabulated. A calculation in any specific 
model boils down to combining the standard loops and using the tables. A 
number of recipes for the reduction to minimal operators can be found in [24j . 
Doing loops becomes a business similar to doing integrals. 

The fact that some coefficients in the effective action remain unspecified 
is none of the tragedy. The effective action is a phenomenological object in- 
tended for obtaining the values of observables. The spectral weights are cer- 
tain phenomenological characteristics of the vacuum like the permittivity of 
a medium. They are to be calculated from a more fundamental microscopic 
theory. Some microscopic theory of some level is incapable of specifying some 
of the coefficients. So what? Classical theory was capable of even less, and, 
nevertheless, celestial mechanics has been successfully worked upH. The only 
important question is whether the lack of knowledge affects the problems that 
we want to solve. This will be cleared up in the next lecture. 



Remarkably, without a knowledge of string theory! 
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4 Vacuum Currents and The Effect of Particle Creation 



Vacuum Currents 

Consider quantum electrodynamics. In this case, (x) is a set of the vector 
connection field and the electron-positron field 

QED: ^'^=(^^,V'). (4.1) 

The commutator curvature is, up to a coefficient, the Maxwell tensor, and the 
operator field equations are of the form 

^''Tl^M) + JM = -JT (4-2) 

where Jju(V') is the operator electron-positron current, and J'^*' is an exter- 
nal source. Averaging these equations over the in-vacuum state, one obtains, 
according to the general derivation above, the same terms but as functions of 
the mean field plus a set of loops: 

^''n.,{{A)) + J,{J^)+^(/^+^^+^(^ (4.3) 



A /A / ip yj 

There is another such equation, for ip, but, since ip has no external source, its 
solution is 

(V) = . (4.4) 

Then, in (4.3), J^d'tp)) vanishes, and the loops with the vertices S^^^ vanish. 
There are no such vertices in QED but, if there were, as in gravidynamics, 
they would be proportional to (tp) and vanish by (4.4). The photon loop also 
vanishes because neither there is a vertex Sj\^_aa but this is already a specific 
property of QED. Only the electron-positron loop survives. 

The surviving loop is a function of (A) , and, by derivation, is the electron- 
positron current averaged over the in-vacuum: 

\j = Jrm) = (in vac| J^(V')|in vac) . (4.5) 

This is the vacuum current. According to (4.3), the observable electromagnetic 
field satisfies the Maxwell equations with an addition of the vacuum current: 

w^n^^iA) = -j;'^iA) - . (4.6) 

We obtain this current by varying the effective action and next replacing the 
Euclidean resolvents with the retarded resolvents: 



T-vac/ ^\ _ Sr{A) 

^ ^ '~ 6A^ 



(4.7) 
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r{A) = J dxg^l^ \nF{a)n + F{a^, Da, □3)7^l7^27^3 + . . .] . (4.8) 

It is completely similar if ^p°^ (x) is a set of the metric field and any matter 
fields 

GRAVITY: = (3^,, V-) • (4.9) 

The only difference is that the vertex Sggg is nonvanishing: 

V' 9 

RA{9)) - 7;{9^.)Ri{9))+H^ + ^{^ = 8nT^T , (4.10) 



9 

(V) = , (4.11) 

and it is assumed again that the matter fields have no sources. Again, by 
derivation, the matter loop is the energy-momentum tensor of the field 
averaged over the in-vacuum but the vacuum current contains, in addition, 
the graviton loop: 

T^^= == (in vac|T^^(V')|in vac) + the graviton loop. (4.12) 

The Einstein equations are replaced by the expectation- value equations in the 
in-vacuum state: 

R^. - \9;..R = 87rT;-(5) + SttT-* . (4.13) 

Since the gravitational field couples to everything, the equation (4.10) 
should contain loops of all matter fields in Nature. The effective actions for 
all loops including the graviton loop have the same structure: 



(4.14) 



r{g) = J dxg'/^ [RFia)R, + F{a,, O^, □3)i?i..i?2..i?3.. + •••]• (4.15) 

Only the coefficients of the formfactors are different. To have the correct 
coefficients, one would need to know the full spectrum of particles. Therefore, 
in the case of gravity, the axiomatic approach is most suitable. 

Now recall that the curvatures are redundant, and the effective action is 
in fact a functional of the conserved currents (3.16) and (3.17). Owing to this 
fact, the expectation-value equations (4.6) and (4.13) close with respect to 
these currents: 

(v^Te,^) + /(□ret)(v''7e,^) + o(v^7^,^)' = -j;^' , (4.i6) 
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+ /2(aret)(V^V. - g^.a)R + o(r^^ - ^g^.nf = SttT^^* . (4.17) 

Of course, with respect to the mean fields, these equations are closed from 

the outset but, at an intermediate stage, they are closed with respect to the 
Maxwell and Einstein currents. When solved with respect to these currents, 
they become literally the Maxwell and Einstein equations with some external 
sources but not the original ones. To make this clear, use the fact that the 
vacuum terms are proportional to the Planck constant and solve the equations 
by iteration: 

v^Te.^ == -r^' + /(Dret) + o {JTf . (4-18) 

R^.^ - Ig^^.R = 8^T°f - ./i(nrct)8^r-^ 

+ ./2(nrct)(V^V, - 5p.n)8^r'=''* + O {T;ff . (4.19) 

These arc the Maxwell and Einstein equations with the original sources prop- 
agated in a nonlocal and nonlinear manner. 

There is an effect in these equations that drives the entire problem. 

Emission of Charges 

Consider again QED and suppose that the external source has a compact spa- 
tial support. This source is the current of a set of electrically charged particles 
moving inside a spacetime tube but, since the observable electromagnetic field 
is the expectation value, only the total current in (4.6) or (4.18) is observable: 

jtot ^ jext + jvac(_4^ _ (4 ^0) 

And the total current has a noncompact spatial support because the vacuum 

contribution is nonlocal. One may calculate the flux of charge through the 
support tube of J*^'^* and even through a wider tube (see Fig. 3), and it will 
be nonvanishing: 

er{Si) - er{S2) = dTi^j^O. (4.21) 

Here eq-i^) is the amount of the electric charge contained inside the tube T 
at a given instant S. The charge inside the tube is not conserved. 

If, when moving away from the support of J^^^, the flux (4.21) falls off 
rapidly, then its nonvanishing only means that the boundary of the original 
source gets spread. Because of the creation of virtual pairs, this boundary can 
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— supp J'^ 




Fig. 3. Support tube of J'^^ and a wider tube. 



never be located precisely. The charges of the external source immersed in 
the quantum vacuum are always annihilated and created again in a slightly 
different place. There is no point to worry about. Just step aside a little. 

However, one may ask if there is a flux of charge through an infinitely wide 
tube: 

eiEi) - e{E2) = ^ ! J""' dT^' ■ (4-22) 



47r 



In this equation, e{S) is the total amount of the electric charge in the compact 
domain of space at a given instant S. For (4.22) to be nonvanishing, 
should behave as 

J,™'^ = O ( 4 1 , r ^ oo , (4.23) 



vac 



V area of S 



(4.24) 



where S is the intersection of T with S (Fig. 3). In this case, it would turn 
out that the charge disappears, i.e., our source is emitting charge. But even 
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this may not be a point of concern if the current in (4.22) oscillates with time, 
and the oscillations sum to zero for a sufficiently long period between Ei and 
X'2. The expectation values have uncertainties, and these oscillations are a 
quantum noise. Just do not measure (4.22) too often. 

However, one may ask if the charge emitted for the entire history 

e(-oo) - e(+oo) = / j;-rfT'* (4.25) 

47r J r— >C!0 

is nonvanishing. There will always be oscillations in the current but they may 
sum not to zero. Since, as r ^ 00, all fields fall off, there are, in this limit, 
the asymptotic Killing vectors corresponding to all the symmetries of flat 
and empty spacctime. Therefore, one may ask the same questions about the 
emission of energy and any other charges. Thus the quantity 

M(-oo) - M(+oo) = / TZ'S,''dT^' (4.26) 

1:^-00 

with the asymptotic timelike Killing vector is the energy emitted by the 
source for the entire history. 

If the total emitted charges are nonvanishing, then this is the real effect, 
and then the question emerges: what are the carriers of these charges? There 
should be some real agents carrying them away. But the particles of the orig- 
inal source stay in the tube. Besides them, there is only the electron-positron 
field but it is in the in-vacuum state. This means that, at least initially, there 
are neither electrons nor positrons. There remains to be assumed a miracle: 
that either the real electrons or the real positrons - depending on the sign of 
the emitted charge - get created. Then they are created by pairs, and, say, the 
created positron is emitted while the created electron stays in the compact 
domain. 

This crazy guess can be checked. We have two ways of calculating the 

vacuum currents: through the effective action and by a direct averaging of 
the operator currents as in (4.5) and (4.12). Specifically, for the in-vacuum of 
electrons and positrons we have 

T;- == (in vac|T^,('^)|in vac) (4.27) 

where Tf^u{t/j) is the operator energy-momentum tensor of the electron- 
positron field ip. The equation for tp 

(p + li-iq{4))^P = Q (4.28) 

contains the clcctroniagncitic field which in (4.27) figures as an external field 
but is in fact the mean field solving the expectation- value equations. We know 
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that, in the past, all mean fields are static. In the future, they become static 
again because, if the total emitted charges are finite, then all the processes 
should die down. Thus, there are two asymptotically static regions: in the past 
and in the future. The carriers of the emitted charges should be detectable in 
the future as particles with definite energies. But then the state in which they 
are absent is the out-vacuum whereas their quantum state is the in-vacuum. 
It may be the case that the in-vacuum contains the out-particles. This will be 
the case if, between the static regions in the past and future, there is a region 
where {A) is nonstatic because then the basis functions of the Fock modes 
that are the eigenfunctions of the energy operator in the future and the basis 
functions that are such in the past are different solutions of the Dirac equation 
(4.28). 

If we expand 'tp in the basis solutions of the out-particles, insert this ex- 
pansion in (4.27), and then insert (4.27) in (4.26), the result will be 

M(-oo) - M(+(X)) = (in vacl^e^ a+^j^aout^ in vac^ (4.29) 

A 

where ea is the energy of the out-mode A, and similarly for the other charges. 
This result needs no comments. Miracles happen. 



Emission of Charges (Continued) 

An important point concerning miracles is that they happen not always. Let 
us see what is needed for this particular miracle to happen. For that, it is 
necessary to introduce characteristic parameters of the problem. There are 

two sets of parameters. 

Parameters of the quantum field: q, /i. 
Parameters oj the external source: e, I, v. 

Here, q and \x are the charge and mass of the vacuum particles (e.g., of the 
electrons and positrons), e is the charge of the external source, I is the char- 
acteristic width of its support tube, and v is the frequency parameter that 
characterizes the nonstationarity of the source. 
The vacuum current in (4.18) is of the form 

oo 

= [ drn^ p{m^)^^ r^"^ -\- O ( J*^^*)^ . (4.30) 



Here and above, the notation nrct is to record that the resolvent is to be taken 
retarded. The structure of the nonlinear terms in (4.30) is similar: there is an 
overall resolvent acting on a function quadratic in J^^^ (see (2.50)). If the 
vacuum particles are massive, the spectral weight will be proportional to the 
^-function: 
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p{m^) cx e{m^ - (4.31) 

to tell us that there is a threshold of pair creation. We need to find the 
behaviour of J™'^ at a large distance from the support of J^'^': 



r 



=? (4.32) 



First we need to calculate the action of the retarded resolvent on a source 
J'^^^ having a compact spatial support. If J°^^ is static, the result is 



□r 



c 

— exp(-TOr) , J'"'* static. (4.33) 



At a large distance from the source, this is the Yukawa potential. Because 
the function (4.33) is static, it does not depend on the spacetime direction in 
which the limit r ^ / is taken. If J"''* is nonstatic, this is no more the case. 
The limit r 3> Hs direction-dependent, and there are directions in which the 
decrease is slower. Namely, in the directions of the outgoing light rays. 



_ c 



exp (^-m^/rU^ , J°^^ nonstatic (4.34) 



where t7 is a function of timely] whose order of magnitude is 

[/ - i . (4.35) 

V 

Expression (4.34) is to be inserted in the spectral integral (4.30), and, 
since the spectrum is cut off from below, we find that the vacuum current is 
suppressed by the factor 

J™^-exp(^-^^ , r»/. (4.36) 

This is what constrains miracles. However, we find also that the suppressing 
factor depends on the frequency of the source and can he removed by raising 
the frequency. The farther from the support of J'"*^', the greater the frequency 
should be for the current to be noticeable. The pair creation starts as soon as 
the energy Tiv exceeds the threshold 

hv > 2(10^ (4.37) 

but, for the source to emit charge, the frequency should be even greater: 

hv > (pc^) (^/) . (4.38) 

^ Of the retarded time since the surfaces E to which the outgoing light rays belong 
are null. 
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This is easy to understand. The particles start being created in the support 
of the source with small momenta and cannot go far away. The extra factor 
{^c/h)l in (4.38) may be interpreted as the number of created particles for 
which there is room in the support of the source. If the creation is more 
violent, the particles get out of the tube. This is the meaning of condition 
(4.38). The mechanism of emission and conservation of charge is illustrated 
in Fig. 4. There are initially the charges of the external source in its support 
tube. They repel the like particles of the created pairs and, when the number 
of the latter exceeds (^c/h)l, push them out of the tube. The unlike particles 
stay in the tube and diminish its charge. 




Fig. 4. Mechanism of emission and conservation of charge. 



Since the cause of the vacuum instability is the nonstationarity of the 
external source, it is interesting to consider the case where the energy hv 
exceeds overwhelmingly all the other energy parameters of the problem. One 
can then study the strong effect of particle production. It is assumed, in 
particular, that ?iv exceeds both the rest energy of the vacuum particle and 
its Coulomb energy in the external field: 

hv > iic^ , (4.39) 

nv»'^ . (4.40) 

In the limit (4.39), the flux of charge at a given distance from the source ceases 
depending on the mass /z, and the vacuum particles can be considered as mass- 
less. Condition (4.40) enables one to get rid of the consideration of the static 
vacuum polarization which is irrelevant to the problem. The approximation 
(4.39) and (4.40) is called high-frequency approximation. 
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The effective action has been calculated above as an expansion in powers 
of the curvature but the conditions of validity of this expansion have not been 
discussed. This lack can now be met. It is the high-frequency approximation in 
which this expansion is valid. Indeed, consider the series (4.8). Every next term 
in this scries contains an extra power of TZ. and, by dimension, its formfactor 
contains an extra power of D"^. The commutator curvature is proportional 
to the charges and to h~^: 

n^^^ (441) 

In the limit r I along the outgoing light rays, the operator □ contains one 
time derivative: 

□ ~ y . (4.42) 

As a result, every next term of the series contains, as compared to the previous 
one, the extra factor 

In addition, the formfactors in (4.8) can be calculated in the massless limit, 

as has been done above. 

However, the inquest of miracles is not yet completed. Assuming that the 
vacuum particles are massless or that the high-frequency regime holds, we 
get rid of the suppressing exponential in (4.36) but we still need to check the 
power of decrease of the current. The power should be the one in (4.23) for 
the emission of charge to occur. We can readily check this since we know the 
behaviour of the resolvent. Expression (4.34) is again to be inserted in the 
spectral integral (4.30) but this time assuming that the spectrum begins with 
zero mass: 

oo 

= / dm^ p{rn'^)—exp (-mVrU] . (4.44) 
r»! J r \ J 



We see that, for the current to decrease as 0(l/r^), the spectral weight should 
have a finite and nonvanishing limit at zero mass: 

p(0) = finite ^ . (4.45) 

For the respective formfactor, this is a condition on its behaviour at small □. 
The behaviour should be 



F{U) 



OO 



a^o 



We arrive at the following consistency condition on the vacuum formfac- 
tors. In the limit where one (any) of the □ arguments is small and the others 
are fixed, the formfactors should not grow faster than ln(— □): 

F{U) = const. In(-n) , (4.47) 
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F{D,,D2,D3) = /(□2,n3)ln(-ni) , (4.48) 

□ i^O 



If they grow faster, the charges cannot be maintained finite, i.e., an isolated 
system cannot exist in such a vacuum. If they grow as ln(— □), the theory of 
isolated systems is consistent but these systems emit charges. If they grow 
slower, the charges are conserved. 

One can check whether the one-loop formfactors satisfy this consistency 
condition. The second-order formfactors (3.90)-(3.94) do. The third-order 
formfactors behave generally as [35j 

F(Di,n2,n3) = /(Da, □3)^+5(02, Da) In(-ni) + ... . (4.49) 

The alarming terms I/O appear only in the arguments acting on the grav- 
itational curvatures. Therefore, they can affect only the vacuum energy- 
momentum tensor, and it has been checked that, in the energy-momentum 
tensor, these terms coming from different formfactors cancel. In the currents, 
the one-loop formfactors satisfy strictly the consistency condition. Since, in 
addition, their asymptotic ln(— □) terms are nonvanishing, the emission of 
charges in the high-frequency regime is real. The only thing that remains to 
be checked is that this emission is not a pure quantum noise. It will be checked 
by a direct calculation. 

Now one can answer also the question about the indefinite local terms 
in the effective action. The coefficients of these terms are the unspecified 
constants in (3.90)-(3.94). In the limit □ ^ 0, the values of these constants 
are immaterial. Only the terms ln(— □), □ ^ of the formfactors work, and, 
therefore, the incompleteness of local quantum field theory does not affect the 
presently considered problem. 

It will be noted that there are now two mechanisms by which an isolated 
system can emit energy. One is purely classical: a nonstationary source can 
emit the electromagnetic or gravitational waves. The other is quantum: im- 
mersed in the vacuum, a nonstationary source can emit also charged particles. 
A high-frequency source will generally emit both. 



Particle Creation by External Fields 

The problem of particle creation by external fields is a part of the expectation- 
value problem. In the context of the foregoing, it can be set as follows. Consider 
the quantum field that satisfies a linear second-order equation 

(ff^'^V^V J + P) = (4.50) 

containing three external fields: the metric, the connection, and the poten- 
tial. The external fields are asymptotically static in the past and future but 
otherwise arbitrary except that their currents 
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Jap = RaP - ^OapR , (4-51) 

J„ = V'^7^„^ , (4.52) 
Q = P + ^Ri (4.53) 

are confined to a spacetime tube. The quantum field is in the in- vacuum state. 
What is the energy of the quanta of the field (j) created by the external fields for 
the entire history? In the high-frequency approximation, we have everything 
to answer this question. 

To formulate the answer, I need some preliminary construction. Every 
current has an associated quantity called its radiation moment. It will now be 
defined. 

Consider a timelike geodesic in the external metric of equation (4.50). It 
enters the domain of nonstationarity of external fields with a definite energy 
and goes out of this domain with a definite energy. Let E be its energy per 
unit rest mass on going out. I am only interested in the geodesies that escape 
to r = oo. They have E > 1, and, instead of E, I shall use the parameter 7 
defined as 

7=^^^^, E>1, 0<7<1. (4.54) 

At r = 00, the geodesic has a certain spatial direction, or, equivalently, it 
comes to a certain point of the celestial 2-sphere. I shall denote this sphere as 
S, its points as 6: 

e = {61,62), 6eS, (4.55) 
and the integral over the unit 2-sphere as 

d''Si6) (■••). (4.56) 

A geodesic with given 7 and 6 will be called 7, 9 -geodesic (see Fig. 5). 

A. ^,9 -geodesic can be emitted from every point of a compact domain. 
Therefore, the 7, 9 -geodesies with the same values of 7 and 6 make a congru- 
ence, and it can be proven that this congruence is hypersurface-orthogonal. 
Let the orthogonal hypersurfaces be 

T^eix) = const. (4.57) 

Since the parameters 7, 9 fix the congruence, they fix also the family of the 
orthogonal hypersurfaces (4.57), and the "const." in (4.57) fixes a member 
of the family. The fiuiction T^e is determined up to a transformation T-^o — > 
/ {Tj$). This arbitrariness will be removed by the normalization condition 

(VT^,)' = -(l-7') (4.58) 

and the condition that the vector VT-yg is past directed. It is a property of 
the geodetic congruences that the norm in (4.58) can be chosen constant. 




Fig. 5. A 7, -geodesic. 



The radiation moment of any scalar current J is the foUowing hypersurface 
integral: 

D=^ J dx g^/'^S {T^e{x) - r) J{x) . (4.59) 

If the current is not a scalar, it should first be parallel transported from the 
integration point to r = oo along the respective 7, 6 -geodesic. Thus if the 
current is a vector, its radiation moment is 

= [ dx g^/'^S (T^e(x) - r) J^(a;)a^"(a;, 00) (4.60) 
Att J 

where ai3°'{x, 00) is the propagator of parallel transport of vectors to infinity 
along the 7, 9 -geodesic emanating from x. The radiation moment D" is then 
a vector at infinity. In the same way, the radiation moment is defined for any 
current. For the three currents (4.51)-(4.53), the radiation moments will be 
denoted respectively as 

Jaf3, Ja, Q > -Da/3, D^, D . (4.61) 

Since the indices of the radiation moments pertain to a point at infinity, their 
contractions like 

DocD" = QocpD'^D^ , etc. (4.62) 

always assume the flat metric ga^ at infinity. All radiation moments are func- 
tions of four parameters: 

D = D{j,e,T). (4.63) 
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In the limit 7 = 1, the 7, 9 -geodesies become null. The orthogonal hyper- 
surfaces (4.57) also become null, and the geodesies themselves become their 
generators. For the radiation moments, this is a regular limit. Nothing special 
happens to them in this limit except that they become very important. The 
radiation moments at 7 = 1 govern the emission of waves in classical theory. 
Thus if Ja in (4.52) is an electric current, then the following expression: 



(a/(-oo) - M{+oo)^ 



electromagnetic waves 



(4.64) 



7=1 



is the energy of the electromagnetic waves emitted by this current for the 
entire history. A similar expression with the tensor current (4.51): 



(^M(-oo) - Mi+oo)^ 



gravitational waves 



dr 



7=1 



(4.65) 



is the energy of the gravitational waves emitted by the current Jap for the 
entire history. 

The radiation moment is a generating function for the multipole moments. 
The multipole expansion is the expansion of at 7 = 0. It makes sense for 
nonrelativistic systems since 7 is proportional to 1 /c. 

Expressions (4.64) and (4.65) are the solutions of the classical radiation 
problem. And here is the solution of the quantum radiation problem [50] : 



(m(-oo) - M(-hoo)) 



(47r)2 



djj 



created particles 
dr I d^S{9)ir 



dT^ 



D 



1 



9al3 



dr 



1 

"3(1^ 

1 - 1 



dr 
dr^ 



D 



Q/3 



dr2 



(4.66) 



This is the energy of the quanta of the field </> created by the external fields 
for the entire history. As compared to the expressions above, there is an extra 
time derivative in the case of the tensor and scalar moments. It accounts for 
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the dimension of the coupUng constant. Also, instead of setting 7 = 1, one 
needs to integrate over 7. Otherwise, the similarity is striking. The quantum 
problem of particle creation becomes almost the same thing as the classical 
problem of emission of waves. 

The presence in (4.66) of an integral over 7 is not just a technical de- 
tail. The radiation moments have both the longitudinal projections, i.e., the 
projections on the direction of the geodesic at infinity and the transverse pro- 
jections. Inspecting the contractions of the moments in (4.64) - (4.66), one can 
see that, at 7 = 1, the longitudinal projections drop out of these contractions. 
In the integral over 7, also the longitudinal projections survive. Owing to 
this fact, spherically symmetric sources cannot emit waves but can produce 
particles from the vacuum. 

Now I can explain why, when expanding the effective action, I stopped at 
the terms cubic in the curvature. In the high-frequency approximation, the 
expansion (3.40) needs to be calculated up to the lowest-order terms that give 
a nonvanishing effect. The terms of first order in the curvature are local and 
give no effect. The terms of second order in the curvature are nonlocal and 
contribute to the energy flux at infinity but it turns out that their contribution 
is a pure quantum noise. The real effect of particle production begins with 
the third order in the curvature. Expression (4.66) results from the triangular 
loop diagrams. 

Since varying the action destroys one curvature, a cubic action generates 
a quadratic current. This gives the radiation energy a chance to be positive 
definite. Expression (4.66) is positive definite indeed: 

fM(-oo) - M(-l-oo)) , , >0. (4.67) 

V / created particles 

In particular, for the matrix contributions, this follows from relations (3.11), 
(3.12) and the positive definiteness of the matrix uiab' 



tr I -^D I > , tr 



9o:fS 



d.oA (±D^ 



dr J Kdr 



< . (4.68) 



The positivity of the gravitational-field contribution can be proven directly. 



The Backreaction Problem 



The energy emitted by an isolated system (in all forms) should be bounded 
both from below and from above: it should be positive and less than the energy 
stored in the initial state 

< (m(-oo) - M(-hoo)) < M(-oo) . (4.69) 

In expression (4.66), the positivity is guaranteed but the energy c;onservation 
is not. The reason is that the setting of the problem with external fields is 
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physically inconsistent. The vacuum current determines the solution of the 
mean-field equations, and the mean field rather than the external field deter- 
mines the vacuum current. If the backreaction of the vacuum is neglected, the 
conservation laws need not be observed. 

One case in which the vacuum backreaction may not be neglected is where 
both mechanisms of the energy emission, classical and quantum, are engaged 
simultaneously. This concerns particularly the vector connection field. In ex- 
pression (4.66), the integral over 7 has a pole (1 — 7)^^ in the term with the 
vector moment. The residue of the integrand in this pole is precisely the quan- 
tity (4.64), i.e., the energy of the outgoing waves of the vector connection field. 
If it is nonvanishing, e.g., if the external source emits both the electromag- 
netic waves and the electrically charged particles, the integral in 7 diverges. 
The result is a disaster: the radiation energy appears to be infinite. In fact it 
should be taken into account that the created charge affects the generation of 
the electromagnetic waves, and the respective changes in the electromagnetic 
field affect the creation of charge. In the self-consistent solution, the disaster 
is removed. 

Another example concerns the metric field when it has an event horizon. 
In this case, the integral in r diverges at the upper limit. By construction, 
T is the time of an external observer. As t ~^ 00, the source moving in the 
tube hits the event horizon. Its proper time does not turn into infinity. The 
integrand in (4.66) is just finite in this limit, and the integral in r diverges 
linearly. This is the Hawking constant flux of radiation from the black hole. If 
its backreaction on the metric is neglected, the total emitted energy is infinite. 

But even when the quantity (4.66) is finite, it depends on the frequency of 
the source. If the source is external, this frequency is a free parameter. The 
energy of created quanta grows with frequency, and, typically, the ratio 



also grows so that, at a sufficiently high frequency, the energy conservation 
law will be violated. The backreaction should take into account that, when the 
source creates real particles, it loses energy and slows down. It then creates 
less particles, and the process dies away. The conservation laws will then be 
restored. 

The backreaction problem has been solved only in a few cases [5T|-|56j. 
The examples for which it has been solved show that the solution can be 
unexpected and interesting. 

References 



M(-oo) - M{+oo) 
M(-oo) 



In v 



(4.70) 



1. Bryce DeWitt: The Global Approach to Quantum Field Theory, vols 1,2 (Oxford 
University Press, Oxford New York 2003) 



Expectation Values and Vacuum Currents 



55 



2. B.S. DeWitt: Dynamical theory of groups and fields. In: Relativity, Groups and 
Topology. 1963 Les Houches Lectures, cd by C. DeWitt, B.S. DeWitt (Gordon 
and Breach, New York 1964) pp 587 820 

3. G. Jona-Lasinio: Nuovo Cimcnto 34, 1790 (1964) 

4. B.S. DeWitt: Phys. Rep. 19, 295 (1975) 

5. E.S. Pradkin, G.A. Vilkovisky: Lett. Nuovo Cimento 19, 47 (1977) 

6. J. Schwinger: Field theory methods in non-field theory contexts. In: Proc. 1960 
Brandeis Summer School (Brandeis University Press, Brandeis 1960) pp 282- 
285 

7. J. Schwinger: J. Math. Phys. 2, 407 (1961) 

8. L.V. Keldysh: Zh. Eksp. Teor. Fiz. 47, 1515 (1964) 

9. Yu.A. Golfand: Yad. Fiz. 8, 600 (1968) 

10. P. Hajicek: Time-loop formalism in quantum field theory. In: Proc. 2nd Marcel 
Grossmann Meeting on General Relativity (Trieste, 1979), ed by R. Ruffini 
(North Holland, Amsterdam 1982) pp 483-491 

11. E.S. Fradkin, D.M. Gitman: Fortschr. der Phys. 29, 381 (1981) 

12. J.L. Buchbinder, E.S. Fradkin, D.M. Gitman: Fortschr. der Phys. 29, 187 (1981) 

13. R.D. Jordan: Phys. Rev. D 33, 44 (1986) 

14. E. Calzetta, B.L. Hu: Phys. Rev. D 35, 495 (1987) 

15. A.O. Barvinsky, G.A. Vilkovisky: Nucl. Phys. B 282, 163 (1987) 

16. R.C. Hwa, V.L. Teplitz: Homology and Feynman Integrals (Benjamin, New 
York Amsterdam 1966) 

17. G.A. Vilkovisky: Class. Quantum Grav. 9, 895 (1992) 

18. J.S. Schwinger: Phys. Rev. 82, 664 (1951) 

19. J.L. Syngc: Relativity: The General Theory (North Holland, Amsterdam 1960) 

20. G. 't Hooft, M. Vehman: Ann. Inst. Henri Poincare XX, 69 (1974) 

21. P.B. Gilkey: J. Diff. Geom. 10, 601 (1975) 

22. L.S. Brown: Phys. Rev. D 15, 1469 (1977) 

23. L.S. Brown, J.P. Cassidy: Phys. Rev. D 15, 2810 (1977) 

24. A.O. Barvinsky, G.A. Vilkovisky: Phys. Rep. 119, 1 (1985) 

25. G.A. Vilkovisky: Heat kernel: rencontre entre physiciens et mathematiciens. 
In: R.C. P. 25, vol 43 (Publication de I'lnstitut de Recherche Mathematique 
Avancee, Strasbourg 1992) pp 203-224 

26. A.M. Polyakov: Phys. Lett. B 103, 207 (1981) 

27. G.A. Vilkovisky: The Gospel according to DeWitt. In: Quantum Theory of 
Gravity, ed by S.M. Christensen (Hilger, Bristol 1984) pp 169-209 

28. A.A. Ostrovsky, G.A. Vilkovisky: J. Math. Phys. 29, 702 (1988) 

29. l.G. Avramidi: Yad. Fiz. 49, 1185 (1989) 

30. A.O. Barvinsky, G.A. Vilkovisky: Nucl. Phys. B 333, 471 (1990) 

31. A.O. Barvinsky, G.A. Vilkovisky: Nucl. Phys. B 333, 512 (1990) 

32. A.O. Barvinsky, Yu.V. Gusev, G.A. Vilkovisky, V.V. Zhytnikov: J. Math. Phys. 
35, 3525 (1994) 

33. A.O. Barvinsky, Yu.V. Gusev, G.A. Vilkovisky, V.V. Zhytnikov: J. Math. Phys. 

35, 3543 (1994) 

34. A.O. Barvinsky, Yu.V. Gusev, G.A. Vilkovisky, V.V. Zhytnikov: Nucl. Phys. B 
439, 561 (1995) 

35. A.O. Barvinsky, Yu.V. Gusev, V.V. Zhytnikov, G.A. Vilkovisky: Class. Quan- 
tum Grav. 12, 2157 (1995) 



56 Vilkovisky 



36. A.O. Barvinsky, Yu.V. Gusev, G.A. Vilkovisky, V.V. Zhytnikov: Covariant per- 
turbation theory (IV). Third order in the curvature. Report, University of Man- 
itoba, Winnipeg (1993) pp 1-192 

37. A.G. Mirzabekian, G.A. Vilkovisky, V.V. Zhytnikov: Phys. Lett. B 369, 215 
(1996) 

38. Y. Nambu: Phys. Rev. 100, 394 (1955) 

39. N. Nakanishi: Prog. Thoor. Phys. 24, 1275 (1960) 

40. N. Nakanishi: Graph Theory and Feynman Integrals (Gordon and Breach, New 
York 1970) 

41. J. Schwinger: Particles, Sources, and Fields, vol 2 (Addison- Wesley, Reading 
1973) 

42. A. A. Grib, S.G. Mamaycv, V.M. Mostcpanenko: Quantum Effects in Intense 
External Fields (Atomizdat, Moscow 1980) 

43. N.D. Birrell, P.C.W. Davies: Quantum Fields in Curved Space (Cambridge Uni- 
versity Press, Cambridge 1982) 

44. N.M.J. Woodhouse: Phys. Rev. Lett. 36, 999 (1976) 

45. A.G. Mirzabekian, G.A. Vilkovisky: Phys. Lett. B 317, 517 (1993) 

46. A.G. Mirzabekian: Zh. Eksp. Teor. Fiz. 106, 5 (1994) [Engl, trans.: JETP 79, 
1 (1994)] 

47. A.G. Mirzabekian, G.A. Vilkovisky: Phys. Rev. Lett. 75, 3974 (1995) 

48. A.G. Mirzabekian, G.A. Vilkovisky: Class. Quantum Grav. 12, 2173 (1995) 

49. A.G. Mirzabekian, G.A. Vilkovisky: Phys. Lett. B 414, 123 (1997) 

50. A.G. Mirzabekian, G.A. Vilkovisky: Ann. Phys. 270, 391 (1998) 

51. G.A. Vilkovisky: Phys. Rev. D 60, 065012 (1999) 

52. G.A. Vilkovisky: Phys. Rev. Lett. 83, 2297 (1999) 

53. R. Pcttorino, G.A. Vilkovisky: Ann. Phys. 292, 107 (2001) 

54. G.A. Vilkovisky: Ann. Phys. 321, 2717 (2006) 

55. G.A. Vilkovisky: Phys. Lett. B 634, 456 (2006) 

56. G.A. Vilkovisky: Phys. Lett. B 638, 523 (2006) 



